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ABSTRACT: We present a calculation of two-point correlation functions of the stress-energy 
tensor in the strongly-coupled, confining gauge theory which is holographically dual to 
the AdS soliton geometry. The fact that the AdS soliton smoothly caps off at a certain 
point along the holographic direction, ensures that these correlators are dominated by 
quasinormal mode contributions and thus show an exponential decay in position space. In 
order to study such a field theory on a curved spacetime, we foliate the six-dimensional 
AdS soliton with a Schwarzschild black hole. Via gauge/gravity duality, this new geometry 
describes a confining field theory with supersymmetry breaking boundary conditions on 
a non-dynamical Schwarzschild black hole background. We also calculate stress-energy 
correlators for this setting, thus demonstrating exponentially damped heat transport. This 
analysis is valid in the confined phase. We model a deconfinement transition by explicitly 
demonstrating a classical instability of Gregory-Laflamme-type of this bulk spacetime. 
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1 Introduction 

The AdS/CFT conjecture [1-3] has become a main tool to study strongly coupled gauge 
theories from the point of view of a dual description in terms of a suitable low energy limit 
of string theory. In particular, the dictionary of gauge/gravity duality provides a way to 
calculate correlation functions of field theory operators from the gravitational dynamics 
of the bulk. The stress-energy tensor of the field theory is induced by the asymptotic 
behavior of the bulk metric itself: The propagation of small perturbations of the quantum 
stress tensor can be studied by solving the problem of graviton propagation in the bulk 
spacetime. See [4] for the calculation of stress-energy tensor correlation functions in strongly 
coupled J\f = 4 supersymmetric Yang-Mills (SYM) theory, and [5, 6] for the generalization 
to thermal N = 4 SYM theory. The poles of these retarded Green's functions are usually 
easier to calculate because they are just the quasinormal modes in the language of bulk 
gravity [6, 7]. 

Besides understanding supersymmetric and conformal field theories (in particular N = 
4 SYM), one would like to describe strongly coupled gauge theories that share more prop- 
erties with (large N) QCD. We will use an approach where we address the issue of finding 
bulk geometries which are appropriate to break supersymmetry and conformal invariance. 
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Supersymmetry may be broken by considering a bulk which contains a Scherk-Schwarz com- 
pactified dimension and imposing antiperiodic boundary conditions on the fermions [8]. It 
has been proposed to consider the AdS soliton as a bulk which is only locally asymptotically 
AdS and has a compact dimension (an S 1 circle) [10, 13]. The introduction of a scale (the 
size of the compact dimension) also allows to break conformal invariance. The AdS soliton 
is constructed by a double analytic continuation of the planar AdS black hole such that 
a compactification takes place as one makes the original time coordinate Euclidean. This 
geometry is horizon-free and it may be visualized as the surface of one half of a cigar which 
caps off smoothly at its tip (infrared floor). The entropy density therefore vanishes in the 
classical limit and the dual field theory is in a confined phase. See [14] for a discussion 
along these lines, and [15] for an approach towards confinement. The AdS soliton is per- 
turbatively stable and has the lowest energy in its class of spacetimes with the particular 
locally asymptotically AdS boundary conditions. 

One expects that the finite size of the geometry along the holographic direction induces 
a quantization of the graviton modes. These quasinormal frequencies appear as the poles of 
the retarded Green's function in the dual quantum theory, as one would expect on general 
grounds [6]; see also [17, 18]. We will explicitly calculate field theory correlators which 
describe momentum diffusion and we will see how the quantization of graviton modes leads 
to an exponential damping of such transport phenomena. 

After understanding some of these properties of the field theory dual to the AdS soliton, 
one can ask what happens if this theory lives on a non-trivial background spacetime. An 
interesting feature of the AdS soliton is that it can be foliated with arbitrary Ricci flat 
slices which take the role of the boundary geometry in holography. As one foliates the 
AdS soliton with a Schwarzschild black hole, one obtains a six-dimensional spacetime with 
two independent physical scales: the Schwarzschild radius r s and the size of the compact 
dimension L T . This Schwarzschild-black string AdS soliton (Schwarzschild soliton string for 
short) serves as a simple model of a bulk geometry that describes a non-supersymmetric, 
strongly coupled plasma around a (non-dynamical) black hole background. The geometric 
parameter r s corresponds to the plasma temperature T oc r s . We will address the question 
of how this field theory propagates thermal excitations near the black hole horizon towards 
infinity. Analogous considerations as in the case of a flat background lead us to expect an 
exponential damping. 

Besides our interest in stress tensor correlators on a black hole background, another 
closely related aspect of these field theories is their deconfinement transition as one lowers T 
[16]. Holographically, this phase transition can be understood in terms of a Hawking-Page 
transition between the Schwarzschild soliton string and the planar AdS black hole as one 
varies r s at fixed L T . In order to determine the transition temperature, we will calculate at 
which point the Schwarzschild soliton string becomes unstable against small perturbations. 
This computation is closely related to the calculation of thermal Green's functions because 
it also involves solving the bulk graviton equations of motion. Qualitatively, one expects to 
find a Gregory-Laflamme-type instability [19] when the black hole horizon is small compared 
to the compact dimension. 

This paper is organized as follows. In section 2, we use the gauge/gravity duality to 
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calculate a stress-energy tensor two-point function in the strongly coupled, confining gauge 
theory which is dual to the AdS soliton. In section 3, we introduce the Schwarzschild soliton 
string and confirm quantitatively that it is indeed unstable against small perturbations with 
tensor modes of a decomposition of the spacetime with respect to the base manifold *B = 
Schw4 (i.e. the Schwarzschild black hole). The instability problem reduces to a combination 
of the well-known Gregory-Laflamme instability and the propagation of a scalar in the 
AdS soliton. Once we have identified the stable phase of the Schwarzschild soliton string, 
section 4 will be concerned with the study of stress-energy correlation functions in the 
strongly coupled field theory on a Schwarzschild black hole background which is dual to 
the Schwarzschild soliton string in gravity. An important ingredient for this calculation 
will be the results from section 2. Due to the fact that the boundary metric is no longer 
translationally invariant, we will need to generalize the Fourier decomposition that could 
otherwise be used. We conclude with some remarks in section 5. In appendix A, details of 
the calculation of the AdS soliton quasinormal modes are outlined. We show in appendix 
B that vector and scalar perturbations do not destabilize the Schwarzschild soliton string. 



2 Shear Diffusion in the AdS Soliton Dual 

In this section we review some properties of the AdS soliton geometry and its dual confin- 
ing field theory. We calculate the quasinormal modes (QNM) for a linearized perturbation 
of this geometry which propagates like a scalar field. The QNM contribution with longest 
wavelength dominates stress-energy correlators in the dual field theory. Besides being inter- 
esting for their own sake, results of this analysis will be needed for calculations in section 4. 
By viewing our quantum theory as a toy model for QCD, the poles of the Green's functions 
may be interpreted as glueball masses. Such holographic computations of glueball spectra 
in QCD3 and QCD4 have been carried out, see e.g. [9]. The authors of [10] even made use 
of the same AdS soliton geometry and calculated glueball masses with a WKB approach. 
We will make a comparison with their results. 

There have also been proposals for models which are more driven by phenomenology. 
In [11] the bulk is constructed as a stack of flat slices which are conformally rescaled in such 
a way as to reproduce known QCD phenomenology. This model relies on similar computa- 
tions involving scalar field propagation and QNMs as the model that we will be studying. 
In this sense our calculation might be easily adopted to such more phenomenologically 
relevant scenarios. Models based on D-branes are also able to reproduce within certain 
bounds the phenomenology of QCD. In particular, the Sakai-Sugimoto ansatz models QCD 
holographically by investigating stacks of D-branes, see [12] and references therein. 
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2.1 The AdS Soliton 

The AdS soliton has been described by Horowitz and Myers in [13]. To construct the 
AdS soliton, we start with the following AdS black hole solution to d-dimensional Einstein 
gravity with cosmological constant A < 0: 



^.2 

ds2 = ¥ 



i - + {dx i ? 



r" 



„d-l 



-1 



where t = — (d — l)(d — 2)/2A is the AdS^ radius, and i = 1, . . . , d — 2. If we now perform 
a double analytic continuation of this metric, i.e. t — > it and x d ~ 2 — > it, we obtain 



ds2 = P 



i hn ,d.v"<lr" - j 1 - ] dr 2 



+ f 1 " 1 (2-2) 



where 77^ is the (d — 2)-dimensional Minkowski metric, and (x^ 1 ) = (i, x 1 , . . . ,x d ~ 3 ). Note 
that the coordinate r has to be periodically identified in a Kaluza-Klein spirit in order to 
avoid a conical singularity at r = ro, i.e. r ~ r + 4ir£ 2 /(d — l)ro. We anticipate already at 
this point the following nice feature of the geometry (2.2): The flat space metric 77^ can 
be replaced by any Ricci flat manifold and (2.2) will still be a solution of Einstein gravity. 
Since this part of the metric corresponds to the non-compact boundary dimensions, we will 
be able to make a transition to boundary theories on a curved background (see section 3). 

Since the r-circle closes smoothly at r = ro, the geometry just ends there and the entire 
spacetime (2.2) is horizon-free and everywhere smooth. At r = ro there is no singularity but 
an infrared floor where the spacetime ends in a cigar shaped geometry. The AdS soliton has 
a translational symmetry along the [d — 2) Minkowski coordinates, and a U(l) symmetry 
along the compact dimension. Due to the periodicity in r the AdS soliton is only locally 
asymptotically AdS. 

What does the dual field theory look like? First of all, the compact dimension allows 
for supersymmetry breaking by means of imposing antiperiodic boundary conditions for the 
fermions along the r-circle. This introduces a mass gap in the fermionic spectrum such that 
the massive excitations decouple in the low energy effective theory. The compact dimension 
also breaks conformal invariance which can be qualitatively understood by the fact that the 
decoupling of massive fermion modes changes the /3-function (see also [15]). Since the AdS 
soliton geometry is horizon-free, the entropy of the AdS soliton vanishes to first order in 
A^ 2 (i.e. in the classical limit), as one would expect for a field theory in a confined phase. A 
first order confinement-deconfinement phase transition at a certain temperature Td ec . > 
is expected to happen in the field theory [14]. On the gravity side, this corresponds to a 
Hawking-Page transition between the AdS soliton and a Schwarzschild-AdS black hole [20]. 

2.2 Energy-Momentum Correlators: Analytic Approach 

We now want to use gauge/gravity duality to calculate energy- momentum two-point func- 
tions in the boundary field theory of the AdS soliton. This analysis can be done in some 
analogy to [21, 22] since the boundary metric is flat and we can therefore use the same 
simplifications that are used to calculate correlators in thermal N = 4 SYM. 
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In case of the AdS soliton we can consider fluctuations of (f> = h\ that propagate in the 
x 3 -direction. This will eventually allow us to calculate the field theory correlator (T^T^) 
holographically. By the same reasoning as in the case of thermal M = 4 SYM [21], the 
remaining 0(2) symmetry of the background metric ensures that (j) decouples and satisfies 
a massless scalar wave equation in the bulk metric. Rescaling the coordinates as z = r^/r 
and y = rQr/£ 2 , the AdS soliton metric (2.2) becomes 



ds 2 



e_ r 

z 2 



a2( ^ s Mink. 



+ (1 - z d - l )dy 2 + (1 - z d - l Y x dz 2 



(2.3) 



where a = r^/i 2 and d^ink is the line element of (d — 2)-dimensional Minkowski space. 



From the massless scalar wave equation in this metric, d a ( 
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where □ is the wave operator on (d — 2) -dimensional Minkowski space. 

From now on, we concentrate on d = 6, although generalizations are straightforward. 
The case d = 6 describes the dual of a gauge theory in 3 + 1 Minkowski spacetime (times 
S l ) and is therefore particularly interesting. We make the additional assumption that the 
solution is independent of y, i.e. homogeneous along the compact circle. This assumption is 
well justified for the long wavelength limit that we are mainly interested in. We now make 
the Fourier ansatz 



{t,x 3 ,z) 



iuit+iqx s 



with 4> q (0) = 1 



(2.5) 



such that (t>(Q\(q) is the Fourier transform of the boundary value (p(t, x 3 ,0), i.e. we demand 
the normalization (f)q{z = 0) = 1. Eq. (2.4) thus yields the mode equation 
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(2.6) 



where (q^) = (w, 0, 0, q) in the zero frequency limit (i.e. u = 0) such that m 2 = —q 2 is the 
boundary mass. 1 This equation has the following fundamental power series solutions near 
z = 0: 
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(2.7) 
(2.8) 



where the recursion relations for the coefficients a n and b n can be found in appendix A. 

It has been argued that in real time thermal AdS/CFT the incoming wave boundary 
condition at the horizon should be used to single out a unique solution [21]. However, in 



1 Note that we will frequently call the eigenmodes of Eq. (2.6) quasinormal modes despite the fact that 
they are actually momentum eigenvalues at zero frequency. 
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the case of the AdS soliton, the solution that we find near z = 1 is not of the form of an 
incoming or outgoing wave. This is related to the fact that the AdS soliton does not have a 
horizon and we have to impose another boundary condition. As pointed out by Witten [8], 
an important condition that should be imposed for any acceptable solution is the Neumann 
condition d(j) q /dp = at the IR floor z = 1. Here, p is the natural coordinate in which the 
"tip" of the metric at z = 1 looks like the origin of polar coordinates 2 . The above condition 
then just expresses the fact that 4> q is smooth at z = 1. In order to impose this boundary 
condition, we look for a power series solution to Eq. (2.6) near the IR-floor. We find the 
Frobenius solution which looks as follows near z = 1: 

2 4 00 

^ 1} (^1) = 1 + ^(1-^ + ^4(1 --) 2 + ...-E c ™( 1 -^' ( 2 - 9 ) 

n=0 

where the c n are also recursively given in appendix A. There, it is also explained that 
the second independent solution near z = 1 cannot satisfy the above described boundary 
conditions due to a divergent term ~ log(l — z). Switching from (z,y) to the coordinates 
(p,<p) which look like usual two-dimensional polar coordinates with origin at z = 1, one 
can easily verify that the solution (2.9) indeed satisfies the above mentioned Neumann 
condition. 

Knowing that (2.9) is a good solution near z = 1 and that any solution near z = 
can be expressed as a linear combination of the solutions (2.7) and (2.8), we need to find 
out what the global solution is. We thus write the solution <p q ^\z) satisfying the Neumann 
condition at z = 1 in the basis of the two fundamental solutions near the boundary: 

$Hz)=A-$d*)+B-<l>%) (2.10) 

with connection coefficients A, B which might depend on q 2 /a 2 but not on z. Near the 
boundary, ^r°] and <j) q a \ have the forms (2.7, 2.8). In Eq. (2.10) we kept the normalization 

(f)q\l) = 1. This could be changed arbitrarily, but as we will see, the Green's function will 
only depend on the ratio B/A, so it would not be affected by another normalization. 

We can now use the same prescription as in the case of thermal M = 4 SYM in order to 
to calculate the stress-energy correlator (T , 2 1 T 2 1 ) [5, 23]. For this purpose we need to write 
down the action of our bulk theory. For approaches to embed the present theory in string 
theory, see e.g. [24-26]. We focus on a universal gauge/gravity duality and consider just 
the low energy gravitational sector which is described by the action 

S = -L f d A xdydz ^g~(K-2A) (2.11) 
J 

with A = -10/ I 2 and the six-dimensional gravitational constant kq. By inserting the 
metric perturbation given by (f>, the part of the (on-shell) action which is quadratic in the 

2 The coordinates (p, tp) in which the IR-floor z — 1 of the d-dimensional AdS soliton (2.3) looks like the 
origin in polar coordinates, are given by p 2 = a(l — z^ 1 ), tp = by with a = A£ 2 /(d — 1) and b = (d — l)/2 
such that tp is 27r-periodic and there is no conical singularity at p = 0. 
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perturbation can be written as 



Squad. = j 4 l o(-q)J r (q,z)(j)o(q)\ Z z= z + contact terms 

1 47T , 

with t = --^—V^gg zz <P- q dz^ q (2.12) 
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with the factor 4tt/5 coming from integrating out the compact dimension. The prescription 
for Lorentzian signature says that we get the retarded Green's function according to the 
following rule [5]: 



G? 2:12 (cj,q) = -2T(q,z) 
_ 27rl 4 a 4 B 
Kit? -A, 



lz->0 

+ contact terms , (2-13) 



where we used the Dirichlet condition <f)q(z = e — > 0) = 1 for the purpose of finding the 
overall normalization. 

The poles of the retarded Green's function are given by the zeros of A. On the other 
hand, setting A = in the matching condition (2.10) would correspond to imposing a 
vanishing Dirichlet condition at the boundary z = 0, which defines just the QNM of the 
AdS soliton geometry. This conforms with the general arguments in [6]. 

If we want to calculate the correlation function in position space, the QNM become 
the essential ingredient because we can replace the Fourier integral by a sum over residues. 
Since we will be mainly interested in the zero frequency limit, we set u = 0, such that 

([T 2 Hx 3 ),T 2 H0)])^,Gf 2il2 (x 3 ) = ?^|^e- 3 ^ g 2 /« 2 )- (2-14) 
Instead of integrating q along the real line, we close the contour with a semicircle in the 

3 

upper complex g-plane. Due to the Fourier exponential e iqx , the arc doesn't contribute to 
the integral and we are left with a sum over QNM residues: 

([T^),Ti(0)]) = -^f>-M* 3 Res 9 4 (2.15) 

° n=l 

The calculation of the QNM and of the residues is decribed in the following paragraphs. As 
mentioned above, the simple poles of 1/^4 are just the QNM. Following the general methods 
in [6], for these particular values of g 2 /a 2 , the expansion (2.9) of is normalizable and 
can be matched smoothly with a linear combination of (ffi \ and <^° 2 over the entire interval 
z G (0, 1); see also [18] for an application of similar methods. The motivation for this is 
the observation that the underlying analytic solution is a power series that converges on 
the entire interval, independent of whether we expand around z = or z = 1. One can 
easily check numerically or by investigating the pole structure of the mode equation (2.6), 
that the radius of convergence of the power series of cf)^ in Eq. (2.9) reaches z = 0, such 
that the connection coefficient A can be found by evaluating the matching equation (2.10) 
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QNM 


Qi 


<?2 


<?3 


<?4 


15 


<?6 


our value/a 


4.061 % 


6.688 i 


9.249 i 


11.786i 


14.313 i 


16.833 i 


WKB estimate/a 


3.96 % 


6.63 i 


9.21 i 


11.75 a 


14.29 i 


16.81 % 


Res 9n B/A 


-21.08 


-244.4 


-1219 


-4072 


-1.07- 10 4 


-2.41 • 10 4 



Table 1. The first two lines show the values of the lowest QNMs in the six-dimensional AdS-soliton, 
using our matching method and the WKB estimate from [10], respectively. The third line shows 
the residues of B/A at these points. 



at z = with the involved functions ($\, </>i°2 and being given by their power series 
expansions: 

oo 

A=J2°n- ( 2 - 16 ) 

n=0 

The discrete set {q n | _4.(g^/a 2 ) = 0} turns out to be purely imaginary. We determine these 
zeros numerically, using partial sums of the explicit expansion of A in Eq. (2.16). A WKB 
estimate for the same eigenvalues has been given in [10], where the set of m 2 = — q^ has 
been associated with the glueball masses in the dual field theory. Their result (rewritten in 
terms of our conventions and parameters) is: 



2 2 

m n = ■ 



i = n(n + iy^(^) +O(n ). (2.17) 



The first QNM are listed in table 1. Even beyond the shown accuracy, the values that 
we obtain from our matching method agree precisely with what we find by just using a finite 
differences algorithm to solve Eq. (2.6) numerically. We observe that the WKB results agree 
with these exact values to an accuracy which is in accordance with Eq. (2.17), becoming 
better for larger n. Note that the complex conjugates of all the q n are also zeros of A. 

The functions <j>q n (z/zo) for n = 3, 5 are plotted in fig. 2.2. The expansion around 
z = is shown on the interval [0,0.95], and the expansion around z = 1 is shown on [0, 1]. 
We cannot distinguish them in the plots since they match perfectly over the entire common 
interval when q G {q n ,Qn}n- For all other values of q the boundary condition (f> q (0) = 
makes it impossible to match the two expansions. 

We can also express Basa function of A. Since all the power series expansions converge 
at z = 1/2, we can evaluate Eq. (2.10) at this point 3 , and find 

B ~ e^) • (2 ' 18) 

We can now calculate the residues of B/A at the QNM. We take Eq. (2.16) and plot 
(q — q n )/A, which is a smooth function in the vicinity of q n , and determine the value at q n 
with high numerical precision. This result is multiplied with the value of B at the particular 
point, B(q^/a 2 ). This latter value can easily be obtained from Eq. (2.18) with the second 



3 Other points in (0, 1) lead to the same result, of course. The convergence might be worse, however. 
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Figure 1. The quasinormal mode functions of the six-dimensional AdS soliton, 4> q , plotted for the 
QNM q$ and q±. The overall scaling has been fixed by normalizing 4> qn {z — 1) = 1 as in Eq. (2.10). 
For the discrete values q € {q n } the expansion around z = and the one around z = z match 
over the entire interval [0, 1). For all other values of q such a matching is not compatible with the 
normalizability condition 4> q (0) = 0. 



summand in the numerator, which is proportional to A, set to zero. This method works 
very well at least for the lower lying QNM. The values of the first six residues are shown in 
the last line of table 1. Using these results, we can evaluate the expression in Eq. (2.15): 



[T}(x%Ti(0)]) 



27^ 4 a 4 



-4.061az 3 



(-21.08) + ... 



(2.19) 



This result confirms our expectation of an exponentially decaying correlation function in 
the long wavelength limit. Physically this means that shear diffusion to infinity is strongly 
supressed. The contribution of the lowest QNM dominates the sum over exponentially 
decaying terms. Although the values of the residues grow (see table 1), the exponentials 
make every higher QNM q n completely insignificant compared to the contribution of q n -\. 
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3 Instability of the Schwarzschild Soliton String 

We will now introduce a modification of the AdS soliton which contains a black hole in 
the boundary metric. Before we calculate Green's functions and holographic transport 
properties in this novel geometry, it will turn out to be useful to carry out a stability 
analysis in terms of linearized perturbations. 



3.1 Generalizations of the AdS Soliton 



It has been noted in [16] that the d-dimensional AdS soliton metric (2.2) can very easily 
be generalized. In fact, one can replace the Minkowski metric r]^ in (2.2) by any Ricci flat 
metric g^i^x^) and the resulting metric is still a solution of Einstein's vacuum equations: 



ds 2 



r 

e- 



g^dx^dx u + 1 




+ 1 




(3.1) 



Again, r ~ r + 4ir£ 2 /(d — l)ro needs to be periodically identified. 

If we choose for g^ v the four-dimensional Schwarzschild metric, Eq. (3.1) describes the 
Schwarzschild soliton string. The geometry looks like Schw4 x S 1 stretched out in a string 
along a AdS radial direction r that caps off smoothly at a finite value r = tq. The two 
relevant physical scales are the Schwarzschild radius r s of the black hole and the radius of 
the compact dimension, i.e. L T = 2£ 2 /5tq. In order to simplify calculations considerably, 
we perform the following transformations: 



z 



e 

ro' 



-y- 



t r s t, 



(3.2) 



where r is the original radial AdS coordinate, and p is the original radial coordinate in 
Schw4. This brings the Schwarzschild soliton string metric in the form 



ds' 



f 2 



a 



1 



dt 2 + I 1 



-i 



df 2 + f 2 dW 



+ (1 - z 5 )dy 2 




where we defined the parameter a = ror s /£ 2 . 

The Schwarzschild soliton string is supposed to describe the bulk dual of a strongly 
coupled field theory in a Schwarzschild black hole background. It inherits the important 
property of the AdS soliton that it caps off smoothly at the IR floor deep in AdS. A more 
naive choice for a gravity dual of a field theory in a black hole background would be the 
AdS black string [28]. However, a serious problem would be that the AdS black string is 
nakedly singular at the end point along the string direction. Also, as discussed in [27], 
trying to cover this singularity by a horizon does not eliminate the presence of nakedly 
singular surfaces. The Schwarzschild soliton string clearly solves these problems: There is 
no naked singularity due to the special geometry inherited from the AdS soliton. 
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3.2 Reduction to Gregory-Laflamme Instability 

We want to show that a decomposition of a perturbation of the Schwarzschild soliton string 
geometry in tensor, vector and scalar modes with respect to the base manifold 23 = Schw4 
produces an instability in the tensor sector. The decomposition of the metric (3.3) reads 

ds 2 = g A Bdx A dx B + a 2 (x A )dsl chw , ds| chw = g^ u dx^dx v , (3.4) 

where ju, v run over the indices of the four-dimensional Schwarzschild metric g^ u in the co- 
ordinates of (3.3), gAB describes the two-dimensional orbit space which is parameterized by 
(y, z), and a(z) = la/z. Let us start with a transverse tracefree (TTF) tensor perturbation 
of the form 

9ab — > 9ab + K b , h ab = jjV h\ = = h a \ b , (3.5) 

where a, b = 0, . . . , 5 run over all coordinates. For a solution of the Einstein equations of 
the form 

ds 2 = a 2 (z)dsl chw + £{z)dy 2 + r](z)dz 2 , (3.6) 
the linearized Einstein equations read 

= {-5 c jiD - 2R a c /) h cd +2R\ a h b)c + 2V {a V c h b)c - V a V b h 



=A L h ab 



a A V 277 

2 

a \a J 2 a \( r\ 



8 h 



+ 

V 



1 -32, 



hpv - ^dyh/jtv , (3.7) 



where Al is the Lichnerowicz Laplacian and Al is the Lichnerowicz operator on the 
four-dimensional base manifold 03. We choose a harmonic dependence on y. Further- 
more, the z-dependence is the same for each component h ab . Writing h^ u {x^ ,y,z) = 
Xf_iv{x ti )e LVy H {z) / z 2 , we find that the z-dependence of this equation separates: 



0= [A L + m 2 ) X ^, (3.8) 



5z 4 4\ „, , „ / m 2 1 u 2 



= «"<*> ~ (1^ + - z ) + (Y^w ) m • (3 ' 9) 

where m 2 is the constant that comes from separating the variable z. For v = (no excitation 
in the compact dimension) the second of these equations is exactly the same as that of a 
scalar propagating in the AdS soliton, Eq. (2.6). We can therefore use the results that we 
derived in section 2.2: There is an infinite tower of discrete values of m 2 jo? > for which 
Eq. (3.9) has a regular solution. The existence of an unstable mode thus depends on the 
existence of such a mode which solves Eq. (3.8). But Eq. (3.8) is just the equation that 
governs perturbations of the five-dimensional Schwarzschild black string, i.e. the well-known 
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Gregory-Laflamme problem [19]. We assume that the time dependence of Xuu is of the form 
e Qt . Using the same numerical methods as for the solution of Eq. (3.9), we can show that 
there exists a threshold mode with f2 = which corresponds to the maximum boundary 
mass ?n max ~ 0.768 such that all < m 2 < m max give modes which grow in time. This is 
in accordance with the results of [19]. 

However, we want to investigate Eq. (3.8) in some more detail since we will need the 
form of the analytic solution in section 4. We use the spherically symmetric threshold 
Gregory-Laflamme ansatz: 



(X 



flUj 



fh (x) hi{x) \ 

hi(x) h 2 (x) 
K(x) 
\ K(x)sm 2 6j 



(3.10) 



where we introduced the coordinate x = 1/f with < x < 1. Imposing the TTF gauge 
condition on this ansatz, the perturbation equations (3.8) reduce to the following set of 
equations which form an effectively one-dimensional problem: 

x 4 (3x 2 - 5x + 2) ti 2 - 2x 3 (3x 2 - 6x + 2) ti 2 - (8x 3 + m 2 (2 - 3x)) h 2 = , (3.11) 

ho(x) = ( 2x ( x ~ W " ( 5x ~ e)^) , (3.12) 

Sx — 2 

K(x) = — *— (-x(x - l)ti 2 + {x- 2)h 2 ) , (3.13) 
6x — 2 

hi(x) = 0. (3.14) 

Solving the first of these equations (e.g. numerically with a finite differences algorithm and 
regular boundary conditions) indeed yields as the only non-negative eigenvalues of this 
system the zero mode m 2 = and the Gregory-Laflamme threshold mode m 2 ~ 0.768. 

3.3 Relation Between Physical Scales 

We can now draw conclusions about the critical values of the size of the compact dimension 
and of the Schwarzschild radius of the boundary black hole. To this end, we need to find 
the simultaneous eigenvalues m 2 of Eq. (3.8) and m 2 /a 2 of Eq. (3.9). Therefore, we need 
to combine the threshold eigenvalue m max ~ 0.768 of the Lichnerowicz operator on the 
background metric with the discrete values for m 2 /a 2 = —q 2 /a 2 from section 2.2. This 
yields 

a= T ^ = J ra™* • g {0.216, 0.131, 0.095, 0.073, . . .} . (3.15) 

" V \i/ a ) 

This list continues and, in fact, it gives an infinite tower of discrete values for a which 
asymptotically approach 0. See fig. 3.3 

The largest value of a corresponds to the most unstable mode. Indeed, since we used 
the critical value m max . for which an instability can occur, the values of a for unstable 
linearized tensor modes cannot exceed the value a C rit. ~ 0.216. Writing this result in terms 
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Figure 2. The values a = ror s /£ 2 for which the perturbation equations (3.8, 3.9) have a solution. 
The leftmost point corresponds to the critical value of a for which an instability sets in. 

of the physically relevant scales r s and the radius of the compact dimension L T = 2£ 2 /5ro, 
we conclude that the critical ratio for an instability to occur is 



It is the interplay between these two parameters which determines the stability of the 
Schwarzschild soliton string. If the horizon radius r s is small compared to the size of the 
compact circle (i.e. not bigger than allowed by the above equation), then the Schwarzschild 
soliton string is unstable. Note that we cannot straightforwardly reverse this reasoning 
without studying the vector and scalar sectors of perturbations. In appendix B we carry out 
the analyses of linearized perturbations in the vector and scalar sectors of the decomposition 
of the Schwarzschild soliton string with respect to *B = Schw4. We demonstrate that neither 
vector nor scalar modes give rise to an instability. 

Holographically, we interpret the instability in terms of the confinement scale of the 
gauge theory. Since we associate the temperature of the field theory with the inverse 
period of the compact direction t, we obtain a temperature for the deconfinement transition 
Tdec. = (l/-^T-)crit. for given r s . This confirms the conjecture of [16] in a quantitative way. 

4 Stress Tensor Correlators from the Schwarzschild Soliton String 

We turn now to the task of calculating stress-energy two-point functions in the boundary 
field theory dual to the Schwarzschild soliton string. The analysis is complicated due to a 




(3.16) 
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lack of translational invariance of the boundary field theory metric. First of all, we need 
again to decide how we use the symmetries of the problem to decompose the perturbation. 
As a starting point, we take the Gregory-Laflamme mode and Eqs. (3.11-3.14) from the 
previous section. This corresponds to a TTF tensor mode with respect to the base manifold 
25 = Schw4, and it will eventually allow us to calculate the correlators 

(TuTu) , (TffTff) , (TggTgg) , (T^T^) , (4.1) 

the first of which is particularly interesting: Since Tu is the energy density, the associated 
two-point function describes heat transport in the field theory. According to the dictio- 
nary of gauge/gravity duality, we need to find the solution to the equation of motion that 
separated into the equations (3.8) (or equivalently Eq. (3.11)) and (3.9). Then we need to 
calculate the on-shell action quadratic in the gravitational perturbation, and take appro- 
priate functional derivatives. 

To our knowledge, an analytic solution to Eq. (3.11) does not exist. However, two 
observations suffice to calculate the desired Green's functions in a suitable limit. First, we 
write the equation in Sturm-Liouville form: 



, ,1-x) 2 A / 8(1 -x) (1 , , , , 

9 * 1 (2 -3s) 2 * 2 * 2 J " U(2-3*) 3 + J(2-3 X y m 1 " 2 = ° • (42 

Then, Sturm-Liouville theory ensures that whatever the exact solution to this equation is, 
in an appropriate normalization the solutions to different eigenvalues m 2 are orthonormal 
with respect to the inner product 

dx^(x)4 mi) 4 m2) = «5 m?mi , (4.3) 

where the weight function w[x) is given by the factor multiplying the eigenvalue in Eq. 
(4.2): 

/ N (1 _ X) , 

i4A) 

This is the first observation. The second observation concerns the asymptotic solution to 
Eq. (3.11). For the solution which is normalizable as x — > and regular near the horizon, 
we find the following asymptotics: 

h W( x ) „ I °° e ~ lmUX x2+H/2 + ° {X)] f ° r X -+ ° ' (4 5) 

2 \d [1- \ (8-m 2 ) (1 - x) + O ((1 - x) 2 )] fori-H, 

We can recover the full tensor mode /i^(x M , y, z) = \^ u {x : 0)e wy H(z) / z 2 by integrating 
over the boundary masses. In particular, the time-independent and y-homogeneous 11- 
component is given by 

d{m 2 ) c (m) h { p (x) , (4.6) 

where H^ m \z) is the function H(z) for a particular eigenvalue to 2 . Note that the corre- 
lators which we will calculate refer to the Schwarzschild soliton string in its stable phase. 
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Therefore the integral (4.6) runs over m 2 > , where m max . is the threshold mode for 

the instability to occur. Eq. (4.6) is the analog of the spatial Fourier transform that we 
used in the translationally invariant case (section 2.2): The functions take the role of 
the orthonormal set of base functions, and the correspond to the boundary values with 
respect to which we will take functional derivatives. We define them in analogy to what 
we did in the case of the AdS soliton, such that the normalization of the "mode" function 
H^ m \z) is given by H^ m \0) = 1. Let us again assume v = 0, i.e. no excitation in the 
compact dimension. Then the Eq. (3.9) for H^ m \z) is the same equation as for the mode 
4>q(z) in the case of the AdS soliton. Thus we conclude from our analysis in section 2.2 that 
the properly normalized solution is 



1 + 2 
6a 2 



+ ^(-m 2 /a 2 )(:"' 



m 

Ua 2 



z 7 + 



(4.7) 



These information can be used to determine the Green's function that is induced by 
h 2 . In principle, we proceed as in the case of thermal M = 4 SYM and the AdS soliton (c.f. 
[21] and section 2.2, respectively). However, the details are more complicated due to a lack 
of translational invariance. Using the equations of motion, the part of the six-dimensional 
Einstein-Hilbert action which is quadratic in h 2 can again be written as a five-dimensional 
integral over the AdS boundary term. We start by identifying the relevant terms in the 
quadratic part of the action (2.11): 

v5 



S, 



quad. 



1 

2k« 



dt dx d6 dcf> dy dz 2£ 4 a 4 sin 



1 



{m'{z) 2 + m(z)H"{z) + 



d{x)h 2 {x) 2 + C 2 {x)h 2 {x)h' 2 {x) + C 3 (x)h' 2 {x) 



(4.8) 



where the dots stand for terms which do not contain the right number of derivatives in z, so 
they will be irrelevant in the eventual application of the gauge/gravity recipe. Furthermore, 
we have defined 



Ci(x) 



2 (9x 2 - 20x + 12) 
x 4 (2 - 3x) 2 



C 2 (x) 



Ax (3x 2 - 7x + 4) 
x 4 (2-3x) 2 



C 3 (x) 



3x 2 (x - if 
4 (2 - 3x) 5 



x 



(4.9) 



We proceed by using partial integration in order to remove all derivatives of h 2 and even- 
tually make the x-integration trivial. The first step consits of integrating by parts the term 
~ (h 2 ) 2 . This yields a new term which contains a second derivative, h 2 . If we express the 
full solution h 2 (x) as an integral over "mode functions" as in Eq. (4.6), the resulting 
integral reads 



'S'quad. 



16vr 2 £ 4 a 4 



dxdzd(m\)d{m 2 2 ) c ( mi ) c (m2) 



d4 mi) 4 m2) + (c 2 - c 3 ) 4 mi) 4 ma) ' 



C 3 h 2 



(mi)^(m 2 )// 



(4.10) 
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where also the trivial integrals over 6, <p, y have been performed, and the time-integral 
has been deleted. The time integral can be omitted because we consider time-independent 
solutions, i.e. the final Green's function will not be localized in time. 

If we replace h^™ 2 ^", using the equation of motion (3.11), there are only terms left which 



look like either ~ h^^h^ 2 ^ or ~ h 2 ' n> h 2 ' 12 " ■ Since the coefficient of the h 2 ' L1> ^'"^""-term 
depends on mi and ?ri2 in the same way as in the definition (4.6), we can absorb the integrals 
over the boundary masses such that this term becomes ~ h 2 h' 2 = d x [(h2) 2 /2] and another 
partial integration in x can be performed. This finally gives an expression that does not 
contain any derivatives of h^™ ■ 



(mi), (m 2 )/ 



(mi), (m 2 )/ 



,(on- S hell) _ 167T 2 £ 4 a 4 



quad. 



5k r 



dxdzd(m\)d{m 2 2 ) c ( m i) c (™2) \-3m 2 2 w{x)h { p\x)h { p\x) 



1 — z*- 



^ff( m i)' ff( m 2)' _|_ 4jy(mi)^j(m 2 )// _|_ . . . 



(4.11) 



The function w(x) is precisely the weight function from Eq. (4.4) for which Sturm-Liouville 
solutions h!^ l \x) are orthonormal. The integral over x thus yields a 5- function in the 
eigenvalue, 5(m 2 — m 2 ). 

We still have to integrate out the z-dependence, and find the boundary value of the 
action. First, we integrate by parts the term ~ ff( m ).£f ( m )", such that it contributes with 
a negative sign to the part ~ (H^') 2 . This produces two more terms, one of which 
is a boundary term that is precisely cancelled by the Gibbons-Hawking contribution to 
the Einstein-Hilbert action (2.11). The second unwanted term is oc H^ m 'H^ m '' and it is 
therefore not relevant for the purpose of applying the gauge/gravity recipe. The remaining 
term ~ (_ff( m )') 2 can be integrated as in the translationally invariant case, and we find 



S, 



(on-shell) 
quad. 



48n 2 ^a 4 

5Kg 



d(m 2 ) (c 



1 



-H i - m \z)H {m ^'{z) 



z=l 



2=0 



contact terms . 

(4.12) 



Following the gauge/gravity recipe, we take functional derivatives and obtain the final result 
for the Green's function in position space : 



G R 



11,11^1,^2 J 



d(m 1 ) d(m 2 ) 
96vr 2 £ 4 a 4 



r2 q (on-shell) 
^quad. ^ 



(mi) %)4 m2) (x 2 ) 

B 



d(m 2 ) m 2 h { 2 m \x 1 )h < 2 m \x 2 ) -A-m 2 /a 2 ) 



A' 

ig2in3 f a4 jr ml h i 2 mn \xi)h i 2 mn \x 2 ) Res, 
n=l 



(-m 2 /a 2 ) , (4.13) 



where we evaluated the integral over m 2 as follows: The poles of B/A lie on the integration 
contour. Therefore, the integral is given by its principal value which can be evaluated 



4 Since we work in Lorentzian signature, we would actually have to follow the recipe [5] and throw away 
the contribution at z = 1. Because this contribution is zero anyway, the result would be the same. 
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by deforming the contour such that it goes around the poles in small semicircles in the 
lower half plane. Closing the contour at infinity, we can apply the Cauchy residue theorem 
similar to what we did in the case of the AdS soliton. In order to determine the asymptotic 
behavior of the bi-tensor G^ bcd (xi, x 2 ), we consider the limiting case that one point lies at 
the horizon, x\ — > 1, and the point where the response is measured is far away, x 2 — > 0. 
For this purpose, we can use the approximate solutions (4.5). From the asymptotic form of 
4 m) (>2 -> 0) in Eq. (4.5), we see that the arc at infinity does not contribute to the Cauchy 
integral over m 2 in Eq. (4.13). We obtain for the retarded Green's function: 

00 , „ 

Gf liU (l, x -> 0+) = N ^2 m 2 n e-l^l/* x 2 +\ m \/ 2 [1 + O(x)} Res mn ( " 

n=l ^ 

^ e -4.06* + ... , (4.14) 

with N = — 192i7r 3 £ 4: a 4: CoCi/ Kq. We can clearly see the exponential decay which means 
that transport to infinity of fluctuations of radial momentum density is strongly supressed. 
Also the dominant role of the QNM {m n } is clearly visible. We have not strictly proven 
the convergence of the sum in Eq. (4.14). But although the residues grow with n, one can 
easily check that for small x the exponential prefactors decay much faster for increasing n. 

We would like to find the same qualitative behavior for heat transport, i.e. for the 
Green's function of energy density correlations, Gq 00 (xi,x 2 ). This can easily be achieved, 
using that ho is given by the action of a first order differential operator D x on h 2 : 

h (x) = D x h 2 {x) =>- G§ )j00 (x 1 ,x 2 ) = D xl D X2 G? hll (x 1 ,x 2 ) , (4.15) 

where the exact form of the operator D x can be read off from Eq. (3.12). We see immediately 
that this kind of transformation preserves the qualitative properties of the Green's function 
and in particular its exponential decay as x 2 — >• + . For definiteness, we nevertheless give 
the result for the other Green's functions in the asymptotic limit: 

G* fl0 (l,x^0 + ) = ivf;|m n | 3 e-l m «l/^ 1+ l m "l/ 2 [l + 0(x)} Res m „ , (4.16) 

n=l ^ ' 

Gf 2)22 (l, x 0+) = N ^ e" |m " l/a; * 1+|m " l/2 [1 + O(x)] Res mn (£) . (4.17) 

n=l 



The exponential decay of correlators of Tq shows that heat transport due to small 
perturbations in energy density near the horizon is exponentially supressed as one goes 
radially towards infinity. 



5 Summary and Discussion 

In order to get a step closer towards studying QCD-like theories via gauge/gravity duality, 
we started by investigating the six-dimensional AdS soliton which is completely smooth 
and horizon-free, but it has one compact dimension which allows to break supersymmetry 
and conformal invariance. In the context of gauge/gravity duality, the AdS soliton serves 
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as a toy model to study a strongly coupled field theory with broken supersymmetry in a 
confined phase. On the other hand, we generalized the AdS soliton by observing that it can 
be foliated along the holographic direction with any Ricci flat metric, in particular with a 
Schwarzschild black hole, giving rise to the Schwarzschild soliton string. The Schwarzschild 
soliton string serves as a toy model to understand the field theory dual to the AdS soliton 
on a non-trivial background spacetime. 

By calculating the correlator (Tg (xi)T2 (a^)) y i a gauge/gravity duality, we have shown 
that momentum and energy diffusion to infinity is exponentially supressed in the quantum 
theory on the AdS soliton boundary. This behavior originates from the fact that the QNM 
appear as the poles of the momentum space Green's function, such that the position space 
Green's functions are dominated by these QNM contributions. The quantization of the 
modes is due to the particular geometry, which caps off smoothly at the IR floor. We leave 
it to future studies to elaborate on other modes of perturbations of the AdS soliton. Qual- 
itatively, one expects very similar results, although the analysis is more involved because 
the vector and scalar perturbations have more than just one non-zero component. 

In the case of the Schwarzschild soliton string we found similar behavior of stress-energy 
correlators. The transport of energy and momentum density from the near horizon region 
towards infinity is exponentially supressed in the boundary field theory. The physical reason 
is again the particular bulk geometry which leads to the dominance of QNM in the Green's 
functions. 

We found that due to the presence of another scale (the Schwarzschild radius r s of 
the black hole), the Schwarzschild soliton string shows an interesting classical behavior 
under small perturbations. The stability analysis of the tensor mode can be reduced to a 
combination of the classical Gregory-Laflamme instability and the propagation of a scalar 
field in the AdS soliton. Depending on the relation between the two physically relevant 
scales (i.e. the size L T of the compact dimension, and r s ), this spacetime is unstable if the 
Schwarzschild radius is small compared to L T . We found this instability at the level of 
spherically symmetric linearized tensor perturbations. We have also shown that the vector 
and scalar modes with respect to the base manifold Schw4 do not develop instabilities. Via 
the holographic duality, this instability is interpreted as a deconfinement transition in the 
field theory. 

It is interesting to explore further holographic duals of field theory states on curved 
spacetimes. In [16], the possibility of black droplets and black funnels has been discussed. 
The black funnel is a solution with a single connected horizon that is dual to the Hartle- 
Hawking state of a strongly coupled plasma around a black hole. Such solutions have 
recently been constructed numerically [31]. Black droplets, on the other hand, have been 
conjectured to describe the final state of the AdS black string instability [32] and are 
still to be constructed explicitly. Understanding the endpoint of the Schwarzschild soliton 
string instability in terms of such solutions might then allow for a complete survey of the 
deconfinement transition in a strongly coupled plasma in terms of different bulk geometries. 
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A Power Series Solution for Scalar Field in the AdS Soliton 

In this appendix, we outline how a power series ansatz leads to the solutions (2.7, 2.8). 

Plugging the ansatz <^q\z) = Y^=oA n z n into Eq. (2.6), we find that the coefficients 
Aq and A§ are free, while all others are given by 



A 1= 0, A 2 = -— 2 A , A 3 = 0, M = —,A 



A 



n+l 



1 

v? — 3n — 4 \a 



,An-\ + (n - 4) 2 ^ n _ 4 for n > 5 



(A.l) 



This yields the two Frobenius solutions (2.7) and (2.8) with exponents and 5, respectively. 
The coefficients a n and b n are defined by the A n by setting either Ao = 1, A$ = or vice 
versa: 



,, in _|^ + ..W( 1+ _|_~= 



z* + ... 



= A Q a nz n I + A 5 ( b n i 



(A.2) 



\n=0 



\n=5 



On the other hand, we can also make an ansatz for a power series solution near the IR 
floor z = 1, i.e. (j)^\z) = "Y^ = QC n (l — z) n . This ansatz gives only one free coefficient Co 
and all others as proportional to it: 



1 

45 



10 



a 



( ', + ( 40 + ^2 ) C 2 



C 4 



1 

80 



IOC) - ( 60 + C 2 + ( 105 + 4 ) C'.i 



a- 



a 



i 



15n 2 - 10n+ ) ( '„ 



20n 2 - 50n + 30 + % I C n _i 

o. 



11+1 ~ 5(n + 1)2 

+ (I5n 2 - 65n + 70) C n _ 2 - (6n 2 - 37n + 57) C n „ 3 + (n - 4) 2 C n _ 4 



(A.3) 



This yields the solution (2.9) near z = 1, where we defined again coefficients that are 
independent of the overall scaling Cq: c n = C n /Co. 
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We find only one solution of this form near z = 1 because the second solution does 
not have the form of a simple power series. The indicial equation at z = 1 has zero as a 
double root, so the second independent solution near the IR- floor contains a term of the 
form ~ log(l — z)(j)q (z). Since this is divergent as z — > 1 and does not satisfy the Neumann 
condition, we discard this solution. 



B Vector and Scalar Perturbations of the Schwarzschild Soliton String 

We want to outline the analysis of linearized vector and scalar channel perturbations of the 
Schwarz-schild soliton string with respect to the base manifold 53 = Schw4. We will show 
that linearized, spherically symmetric vector and scalar perturbations do not give rise to 
an instability. 

B.l Vector Perturbations 

Following the general methods for linearized gravitational perturbation theory as developed 
in [29, 30], we start with a vector harmonic on Schw4 which we write as 

V M = [V (f, 9), Vi(f, 0), 0, 0] , (A + k 2 v )Y^ = = V^V M , (B.l) 

where Greek indices and quantities with a hat refer to 23 = Schw4. Keeping the spherical 
symmetry, this vector harmonic does not have components on the sphere, and its remaining 
components are independent of 4>. Also, the mode is assumed to be time-independent, which 
amounts to considering the threshold mode for which an instability could occur. From 
solving the condition V^V^ = 0, one finds immediately that V\{r, 9) = vi(9)f~ 2 (l — 1/r) . 
We see that this solution is not regular at f = 1. Indeed, also from looking at the defining 
equation for a harmonic vector on 23, (A + ky)Y^ = 0, one can straightforwardly infer 
that for ky ^ 0, V\ must vanish identically. The only remaining non-trivial condition from 
(A + fcy)V M = reads 

^ (cote d e V + d$V ) + ( x " J) (?^V + d 2 f V^j = -k v V , (B.2) 

which is just the eigenvalue equation for the Laplacian on S 2 . Considering one particular 
mode with "angular momentum" I, we write Vo(f, 9) = Pi(cos9) vo(f), and obtain 

ky- l ^)v o+ (l-l)( 2 -v' o+ v'A=0. (B.3) 



The perturbation of the Schwarzschild metric which can be built from the harmonic 
h^ = 2a 2 H T Y^, h Afl = af A V^, h AB = , (B.4) 



vector Vj can be written as: 



with the tensor 



= -^(V^V, + V^) . (B.5) 
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It turns out to be useful to formulate the problem in terms of the variables which are intro- 
duced in [29] for the gauge invariant master equation formalism. A gauge transformation 
adapted to the symmetry of the vector mode is described by a gauge vector £ a = (£„ , £4) 
with = aL(x )V £4 = 0. A gauge invariant variable is then given by 

F A = fA + J^VaHt , (B.6) 

where is the covariant derivative with respect to the two-dimensional orbit space. Fixing 
the gauge with 

£ 2 a 2 

^ = -^H T [V (f,9),0,0,0} , (B.7) 

the variable is effectively set to zero and the perturbation is now TTF. According to Eq. 
(3.7), the linearized Einstein equations reduce for a TTF perturbation to A^^j = 0. From 
this, one can immediately see that F y = such that the Lichnerowicz equation reduces to 
a radial equation in f, and an equation for F z {z). The radial equation is exactly the same 
as (B.3), and the equation for z reads 



±F z (z) + - z F' z (z)-F'J( 



k 2 v 1 



a 2 (1 — z 5 

This equation with regularity boundary conditions allows for discrete modes with ky/a 2 < 
0, very similar to what we found for tensor modes. Since a 2 can by definition not be 
negative, the question is: Does Eq. (B.3) have regular solutions with negative eigenvalues 
ky < 0? If so, then they would correspond to unstable vector modes. However, applying a 
finite differences algorithm, we find that independent of what value I 6 {0, 1,2,.. .} takes, 

„2 
"V 



the equation does not have any eigenvalues ky < that would correspond to an instability. 



B.2 Scalar Perturbations 

The scalar harmonics § are defined by (A + fc|)S = 0. Assuming that our perturbation does 
not break the symmetry of the underlying 2-sphere, i.e. S is an eigenmode of the spherical 
Laplacian and can thus be written as § = Pi(cos6)s(r), the defining equation for S becomes 

We can again find the allowed eigenvalues kg numerically. As it was the case for the vector 
modes, we only find eigenvalues kg > 0. A comparison with vector modes thus already 
stronlgy suggests that there is no instability. However, for definiteness, we complete the 
perturbative analysis in the following paragraphs. 

We will show that all eigenvalues of the perturbation equations satisfiy k s /a 2 < 
which is not compatible with the fact that Eq. (B.9) did not yield any negative eigenvalues 
k s < 0. Therefore the only consistent eigenvalue of the full problem will be kg = such 
that there is no threshold mode and no instability. 

As in [29], we construct the following vectors and tensors build out of the scalar har- 
monic §: 

S M = -^V M S, S^ = ^V^S + ^S, (B.10) 
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where g is the four-dimensional Schwarzschild metric. From these building blocks we can 
construct a symmetry adapted scalar perturbation: 

hfj, v = 2a 2 (H L (z) g^E + H T (z)E f _ iU ) , h Afl = af A (z)S fl , h AB = /ab(z)S . (B.ll) 

We proceed as in the case of vector harmonics. First, we rewrite the degrees of freedom of 
the perturbation h a b in terms of the following (gauge invariant) variables: 

F:=H L + - i —H T + -X7 A fX A , (B.12) 
a — 2 a 

Fab :=fAB + V A X B + V B X A , (B.13) 

where X A := -j^- (f A + ^aH t J • (B.14) 

We fix the gauge by using the symmetry adapted gauge vector 

fa 2 
k s z 



Z = -7TT2 I A i(^, A 2 (z)S, A 3 (z)S] , (B.15) 



with A 1 (z) = H T (z), A 2 (z) = ^X y (z), A 3 (z) = ^X z (z). (B.16) 



This yields the following simple expression for the scalar perturbation: 
h a b 



/V o o i 2(V , 



F m SF yz S\ with h IMU = F(z)—^-Sg^. (B.17) 
y F yz SF zz §; 

We use this perturbation and study the resulting Einstein field equations (3.7). This yields 
a set of equations for F AB {z) and F(z). After some algebraic operations, one finds that 
these equations can be decoupled: 

kln K 2 -2(l + 4z 5 )^F(z)+^z(4 + z 5 )^- 8 ^ ^ 2 

- z 2 (1 - z 5 ) - A (i _ 2 5) ) F "( z ) = o , (B.18) 

U2 



a* \ a 2 



2 (1 + 4z 5 ) W) + (4 + z 5 ) ^ - * (1 - z 5 ) 2 ) F'(z) 

)) *"( 

^) = -2' 2 (l- 5 )(^) + 2(|) 1 (1 - z") d z Q^W)) ' (B.19) 

F„(*) = - 2Z2 (T^5) (^(*) - 2 (3) 1 (1 - ^ (j^)) > ( B - 2 °) 
F yz (z) = 0. (B.21) 

Solving these equations is slightly more complicated than in the vector and tensor sectors 
due to the fact that the eigenvalue k^/a 2 appears quadratically in Eq. (B.18). We solve Eq. 
(B.18) directly by "shooting" the eigenvalue such that the regularity boundary conditions 
are satisfied. We find that (just like in the case of vector and tensor perturbations) there 
are no positive eigenvalues k^/a 2 > 0. Instead, there are infinitely many discrete negative 
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eigenvalues in addition to the eigenvalue 0. Due to the lack of a compatible negative 
eigenvalue of Eq. (B.9), we conclude that at the level of linearized perturbations which 
have the previously assumed form in the scalar sector with respect to Schw4, there are no 
unstable modes. 
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